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^ ; ABSTRACT 

CL( Aims. We present synthetic bolometric and broad-band UBVRI light curves of SNe la, for four selected 3-D deflagration models of thermonu- 
X clear supernovae. 

» I Methods. The light curves are computed with the 1-D hydro code stella, which models (multi-group time-dependent) non-equilibrium radia- 
tive transfer inside SN ejecta. Angle-averaged results from 3-D hydrodynamical explosion simulations with the composition determined in a 
nucleosynthetic postprocessing step served as the input to the radiative transfer model. 

IL" Results. The predicted model UBV light curves do agree reasonably well with the observed ones for SNe la in the range of low to normal 
. 5^ luminosities, although the underlying hydrodynamical explosion models produced only a modest amount of radioactive ^^Ni (i.e. ~ 0.24 - 
0.42 Mq) and relatively low kinetic energy in the explosion (less than 0.7 x 10^' erg). The evolution of predicted B and V fluxes in the model 

5_j with a ^*Ni mass of 0.42 Mq follows the observed decline rate after the maximum very well, although the behavior of fluxes in other filters 
somewhat deviates from observations, and the bolometric decline rate is a bit slow. The material velocity at the photospheric level is of the 
order of lO"* km s"' for all models. Using our models, we check the validity of Amett's rule, relating the peak luminosity to the power of the 
deposited radioactive heating, and we also check the accuracy of the procedure for extracting the '*Ni mass from the observed light curves. 
Conclusions. We find that the comparison between theoretical light curves and observations provides a useful tool to validate SN la models. 
The steps necessary to improve the agreement between theory and observations are set out. 

Key words. Stars: supernovae: general - Hydrodynamics - Radiative transfer - Methods: numerical 



1. Introduction 

With three-dimensional simulations (e.g. Hillebrandt et al., 
2000; Reinecke et al., 2002; Gamezo et al., 2003; Ropke & 
Hillebrandt, 2005b), type la supernova (hereafter SN la) explo- 
sions can be modeled self-consistently. Such models are con- 
structed from first principles, avoiding free parameters in the 
description of physical processes. The only parameters enter- 
ing are the initial state of the exploding white dwarf star and 
the configuration of flame ignition, which must be determined 
in separate studies of progenitor evolution. 

Consequently, the important question is to what degree do 
these models describe real SNe la explosions. This can only 



be answered through comparison with observations of nearby 
events (e.g. Pignata et al. 2004; Benetti et al. 2004; Stehle et al. 
2005). To this end, observables need to be derived from theo- 
retical explosion models. In this paper, we show how synthetic 
light curves (both bolometric and broad-band UBVRI) can be 
obtained from a set of three-dimensional deflagration models 
of SN la explosions (Travaglio et al. 2004; Ropke et al. 2005; 
Ropke & Hillebrandt 2005b). We emphasize that the present 
study is intended to demonstrate the derivation of light curves 
from theoretical models and to analyze the use of these as a tool 
to validate such models. Therefore the chosen models do not 
reflect particularly realistic simulations and have known weak- 
nesses. Judging the validity of explosion models in terms of 
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derived observables requires more elaborate models and this 
issue will be addressed in a forthcoming publication. 

When one succeeds in constructing SN la explosion models 
that are consistent with observations, they will be used to tackle 
questions concerning the application of SNe la in determining 
cosmological parameters. Since these objects need to be cal- 
ibrated on the basis of empirical recipes to serve as distance 
indicators, such distance measurements are subject to uncer- 
tainty. This may be lessened by a theoretical understanding of 
the origin of the diversity in observables and of the caUbration 
techniques. 

The cahbration methods currently applied in SN la cos- 
mology are based on an empirical correlation between the 
blue-band magnitude of the supernova and the decline rate of 
its hght curve (e.g. Pskovskii, 1977; Philhps, 1993; Phillips 
et al., 1999). As pointed out by many authors (e.g. Hoflich 
et al., 1998; Sorokina et al., 2000; Mannucci et al., 2005), the 
luminosity-decline rate relation for SNe la is derived at low 
redshifts z, and it is unclear how it should change at high red- 
shifts. Systematic differences of the order of ~0.1 magnitudes 
(Hoflich et al., 2000; Nomoto et al., 2003) can be important in 
cosmological applications. A sound understanding of the phys- 
ical mechanism of SNe la is crucial in any project investigating 
possible systematic trends in their properties. Additionally, a 
radiative transfer algorithm is required to reliably relate the hy- 
drodynamical model of the explosion to observed fluxes and 
spectra of SNe la. 

Ropke & HiUebrandt (2004) and Ropke et al. (2005) have 
presented systematic studies of eff'ects resulting from variations 
in the progenitor white dwarf's carbon-to-oxygen ratio, its cen- 
tral density at ignition and its metalUcity on 3-D deflagration 
explosion models. These were based on simple and thus only 
weakly exploding models and therefore such studies need to be 
extended with more elaborate simulations. The missing Unk to 
SN la cosmology in such approaches is the derivation of syn- 
thetic hght curves. 

This paper presents synthetic light curves for two models 
from the set described by Ropke et al. (2005) and for two addi- 
tional 3-D models (from Travaglio et al., 2004; Ropke, 2005). 
We compute bolometric and broad-band UBVRI Ught curves 
of SNe la, employing the radiation hydro code stella which 
simulates multi-group time-dependent non-equilibrium radia- 
tive transfer inside SN ejecta in a spherically- symmetric ap- 
proximation. 

First, we briefly outhne the assumptions made in modeling 
the radiative transfer (Sec. 2). In Sec. 3 we describe parame- 
ters of the hydrodynamical models. Various observables pre- 
dicted by the radiation code for the studied hydro models are 
presented in Sec. 4. The accuracy of determining the ^*Ni mass 
based on Arnett's rule is demonstrated in Sec. 5. 

The results are summarized and discussed in Sec. 6. We 
describe the necessary steps for future improvements of both 
hydrodynamical and radiative transfer models in Sec. 7. 

2. Numerical Model for Radiative Transfer 

Modeling the post-explosion stages of SNe la appears to be 
easier than that of core collapse supemovae (CCSNe). For ther- 



monuclear supemovae the hydrodynamics is very simple: the 
coasting stage starts very early, there are no shocks, and the 
only additional heating mechanism is the decay of radioactive 
nuclei. A typical assumption for SNe la Ught curve modeling 
is to fully neglect the hydrodynamics. We discuss the validity 
of this assumption in Sect. 6. 

However, with SNe la difliculties arise with the modeling 
of the radiative transfer. SNe la become almost transparent in 
the continuum at the age of a few weeks. This means that NLTE 
effects are stronger compared to CCSNe. Radiation decouples 
from matter within the entire SN la ejecta prior to maximum 
hght (Eastman 1997; Blinnikov & Sorokina 2004). In this en- 
vironment one cannot ascribe the gas temperature (nor any 
other temperature) to radiation since a SN la spectrum differs 
strongly from a blackbody. Instead, one has to solve a system of 
time-dependent transfer equations in many energy groups with 
an accurate prescription for the treatment of a enormous num- 
ber of spectral hues. These hues are the main source of opacity 
in the ejecta (Baron et al. 1996a; Pinto & Eastman 2000b). 

Recently, powerful codes have been developed in order to 
address the full 3-D time-dependent problem of SN la radia- 
tive transfer (Hoflich 2002; Lucy 2005). Yet there are some ba- 
sic questions, like averaging the line opacity in an expanding 
media, that remain controversial. 

We compute broad-band UBVRI and bolometric light 
curves of SNe la with the multi-energy radiation hydro code 
STELLA (Bhnnikov et al. 1998; Blinnikov & Sorokina 2000). 
The 3-D hydrodynamic models are angle-averaged and used as 
input for the radiation hydro code. Time-dependent equations 
for the angular moments of intensity in fixed frequency bins 
are coupled to the Lagrangian hydro equations and solved im- 
plicitly. Therefore there is no need to ascribe any temperature 
to the radiation; the photon energy distribution may be quite 
arbitrary. 

There are millions of spectral lines that lead to the forma- 
tion of a SN la spectrum, and it is not a trivial problem to find 
a convenient way to treat them even in the static case. The ex- 
pansion makes the problem much more difficult to solve as hun- 
dreds to even thousands of lines give their input into emission 
and absorption at each frequency. 

The effect of line opacity is treated in the current work as an 
expansion opacity according to Eastman & Pinto (1993). The 
Une list is Umited to ~160 thousand entries, selected from the 
strongest down to the weakest lines until the saturation in the 
expression for expansion flux opacity is achieved. The effect of 
extending the line Ust deserves further investigation, and work 
in this direction is underway. 

The ionization and atomic level populations are described 
by Saha-Boltzmann expressions. However, the source function 
is not in complete LTE. 

The source function at wavelength A is 

5., = ethB., + (l-eth)/^, (1) 

where eth is the thermalization parameter, is the Planck func- 
tion, and is the angle mean intensity. In a pure LTE approx- 
imation eth = 1, and in a "pure scattering" treatment eth = 0. 
Baron et al. (1996b) compared the results of a full NLTE treat- 
ment, a pure LTE treatment, and a pure scattering treatment of 
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the lines (in their case described the source function only 
in a line). They found that the pure LTE case (for lines) re- 
produces the overall spectral shape of the NLTE case rather 
well, while the pure scattering case does not. In the NLTE case, 
the collisions within multiplets create a pseudo-thermal pool of 
photons which are in equilibrium among themselves. 

In our code, we have rather wide energy bins that contain 
hundreds of strong lines, so our 5 ^ is a complicated average of 
hues and continuum for a given bin. We can ascribe an arbi- 
trary degree of thermalization to the lines . Even if we treat all 
lines as pure absorbers we have an appreciable contribution of 
photon scattering off electrons, thus the average eth 5^ 1 in a 
bin and 5^ is not in full LTE. Therefore, we prefer to call the 
approximation of full thermalization in lines as 'fully absorp- 
tive lines' . The approximation of the fully absorptive opacity 
in spectral lines allows us to simulate NLTE effects in a simple 
manner. NLTE results (Baron et al. 1996b) and the ETLA ap- 
proach (Pinto & Eastman 2000b) demonstrate that fully absorp- 
tive lines give an acceptable description of the overall spectral 
shape and subsequent optical light curves. 

A similar approach to light curve modehng with LTE opac- 
ities and a comparable number of photon energy groups is 
used by Hoflich et al. (1998) who employed NLTE calculations 
(Hoflich, 1995) to calibrate scattering in hues. A full NLTE 
time-dependent treatment is within reach and is already being 
implemented by other groups (Hoflich 2002; Lucy 2005), but 
much work remains to be done before all problems are solved. 

The deposition of gamma rays produced in radioactive de- 
cays of Fe-group isotopes (mostly due to ^^Ni and ^^Co) should 
be properly considered. After being emitted, the gamma rays 
travel through the ejecta and can end up either thermaUzed or 
in non-coherent scattering process. To determine which way 
they end up one has to solve the transfer equation for gamma 
rays together with hydrodynamical equations. The full system 
of equations should also involve the radiative transfer equations 
in the range from soft X rays to infrared wavelengths for the ex- 
panding medium. 

We do not write down the full set of equations solved by 
STELLA (Blinnikov et al. 1998), but we point out that hydro- 
dynamics coupled to radiation is fully computed (homologous 
expansion is not assumed). Since there are several different ap- 
proaches found in literature, here (and in Sect. 6) we discuss in 
some detail the use of the matter temperature equation. 

The equation for the material temperature T is derived from 
the first law of thermodynamics, 

de ^pdV = de-^dp = dT - I 1^ dp, (2) 

which takes the form: 

oo 

lde\dT , r , , d{r^u)Jdp\ 

i^i ^ = ^ ^ J ^"^'^ - ^^^'^ - '^^1^1 • 

Here p - pip, T) is the material pressure, p is the density, m is 
the Lagrangian coordinate (i.e. the mass inside radius r), Jy is 
equal to ^ J ^ dp.Iy, e is the specific internal energy of matter, 
ay is the absorption coefficient per unit mass, and e is the spe- 
cific power of the local heating (here due to gamma-ray energy 



deposition). In the case of LTE the emission coefficient rjy is 

equal to UyBy. 

We should note that the last equaUty in Eq. (2) is based on 
the expression de + pdV = Tds, where s is specific entropy, 
i.e. not only the first, but also the second law of thermodynam- 
ics for reversible processes is used in the derivation. In other 
words, this relation is vaUd when e and p are functions of p 
and T. While radiation is non-equilibrium in our approxima- 
tion, ionization and atomic level populations are described by 
Saha-Boltzmann equations. Therefore Eqs. (2) and (3) are ap- 
plicable. If ionization and excitation depend on other param- 
eters in NLTE, and/or are described by kinetic equations, one 
must use other forms of the energy equation (see e.g. Sorokina 
et al. 2004). 

Another possible approach is to extract matter energy e 
from the equations of evolution for radiation energy E and full 
energy e + E (e.g. Hoflich et al. 1993). In both cases the matter 
energy e is found with lower accuracy than the radiation energy 
E. This is due to the subtraction of nearly equal terms; either of 
e + E and E in the latter approach or in the heating and cooling 
terms of Eq. (3). 

To find the radioactive energy deposition e we treat the 
gamma-ray opacity as a pure absorptive one, and solve the 
gamma-ray transfer equation in a one-group approximation fol- 
lowing Swartz et al. (1995). The effective opacity assumed is 
Ky - 0.0576 cm^/g, where is the total electron number den- 
sity over baryon density. We have checked that this effective 
opacity gives good agreement for the values of the deposition 
and broad band fluxes found by Monte-Carlo codes after the 
maximum light. The rise part of the light curve is not sensitive 
to the variation of Kj when it is changed by a factor of 3. For 
one of the models (h30, described below) the deposition was 
computed in full 3D-transport by a modification of shdom code 
(Evans, 1998). The agreement with ID version is within a few 
percent. 

The heating by the decays ^^Ni — > ^^Co— > ^^Fe is taken 
into account. It is assumed that positrons, bom in the decays, 
are trapped and deposit their kinetic energy locally, cf. Colgate 
et al. (1980a,b); Chang & Lingenfelter (1993); Ruiz-Lapuente 
(1997); Milne et al. (1997); Ruiz-Lapuente & Spruit (1998) and 
JefFery (1999) on various approaches to this problem important 
for late stages (the 'tail') of SNe la. 

To calculate SNe la light curves our method can use up 
to 200 frequency bins and up to ~400 zones in mass as a 
Lagrangian coordinate on a modest processor. 

Major recent improvements in the code stella (introduced 
after the paper Blinnikov et al., 1998, was published) are de- 
scribed in the following. The ionization of the 15 most abun- 
dant elements (from H up to Ni) is now computed for any stage, 
without introducing an averaged ion approximation. As previ- 
ously mentioned the line list has been extended. When comput- 
ing the opacity we now take into account the change in the com- 
position of the iron peak elements in the decay chain ^^Ni — > 
^^Co— > ^^Fe. The most important advancement is also related 
to the computation of the expansion opacity: instead of interpo- 
lating the opacity in tables for preselected instants of time, it is 
now computed for each time step in each mesh zone. Although 
this approach requires much more processor time than using 
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the tables, it is more reliable when one is interested in finding 
fine details for models with small difl'erences in parameters. 

3. Initial models 

For this study light curves are derived using four hydrodynam- 
ical explosion models. Each model possesses different charac- 
teristics that result in a variation in the computed Ught curves. 

The explosion models are based on a pure deflagration 
scenario. In this case a sub-sonic thermonuclear ("deflagra- 
tion") flame is ignited near the center of a Chandrasekhar mass 
carbon-oxygen white dwarf. Subsequently the flame propa- 
gates outwards generating turbulence due to generic instabil- 
ities. The interaction with turbulent motions accelerates the 
flame propagation velocity such that nuclear burning releases 
a sufficient amount of energy that completely disrupts the 
white dwarf. For a review on SN la models see Hillebrandt & 
Niemeyer (2000). 

Details concerning the numerical techniques used to imple- 
ment this scenario were presented by Niemeyer & Hillebrandt 
(1995), Reinecke et al. (1999a, 2002a) and Ropke (2005). The 
main challenge is to self-consistently model the turbulent com- 
bustion which involves a wide range of scales. To this end a 
Large Eddy Simulation (LES) approach is applied, which di- 
rectly resolves only the largest scales of the problem (i.e. scales 
from the radius of the white dwarf, ~2000 km, down to sev- 
eral kilometers). To account for the effects of turbulence on un- 
resolved scales, a sub-grid scale model is adopted (Niemeyer 
& Hillebrandt 1995). The flame itself is modeled as a sharp 
discontinuity separating the fuel from the ashes applying the 
level-set method (Reinecke et al., 1999a). Since such an ap- 
proach does not resolve the internal flame structure, propaga- 
tion due to burning has to be prescribed. This approach how- 
ever does not lead to any free parameters. The theory of tur- 
bulent combustion states that in the case of burning in the so- 
called flamelet regime, which applies to most parts of a SN la 
explosion, the flame propagation completely decouples from 
the micro-physics. It is solely determined by the turbulent mo- 
tions which can be derived from the adopted sub-grid scale 
model. We emphasize that apart from the initial conditions that 
are determined from the progenitor evolution, this model con- 
tains no free parameters. 

Consequently, the models used to derive the light curves 
from viiry only in the initial conditions, i.e. the composition, 
the central density of the white dwarf, and the manner in which 
the thermonuclear flame is ignited. All models were simulated 
on one spatial octant of the white dwarf and mirror symme- 
try was assumed for the other octants. Ropke & Hillebrandt 
(2005b) showed that this artificial symmetry constraint does 
not obscure the explosion mechanism^ In three of the models, 
the flame was ignited centrally and perturbed from spherical 
symmetry by toroids. Contrary to this, the b30_3d_768 model 
of Travaglio et al. (2004) (in the following denoted as b30), as- 
sumed a multi-spot ignition scenario in which the flame was ig- 

' Asymmetric explosions with possible implications on light curve 
modeling may, however, develop from asymmetries in the flame igni- 
tion configuration. 



Table 1. Model parameters. 



model 


Pc9 


X('*0) 




xC-Ne) 


ZIZq] 


1_3_3 


1.0 


0.38 


0.545 


0.075 


3.0 


2_2_2 


2.6 


0.54 


0.435 


0.025 


1.0 


c3_3d 


2.9 


0.5 


0.475 


0.025 


1.0 


b30 


2.9 


0.5 


0.475 


0.025 


1.0 



nited in 30 small kernels per octant. Evidently, such multi-spot 
scenarios can lead to more vigorous explosions than centrally 
ignited models (Travaglio et al. 2004; Ropke et al. 2005). 

The initial parameters of the models are summarized in 
Table 1. We assumed increasing central densities in the se- 
quence of models 1_3_3, 2_2_2 and c3_3d. Explosion models 
c3_3d and b30 were set up for a 50% carbon-oxygen mixture, 
while the carbon mass in models 2_2_2 and 1_3_3 was set to 
0.46 and 0.62, respectively. 

The detailed nucleosynthesis of each model was derived 
with a post-processing technique that uses data obtained from 
tracer particles in the explosion simulation (see Travaglio et al., 
2004; Ropke et al., 2005, for a detailed discussion concerning 
this method). During the post-processing solar metalUcity was 
assumed for all the models except 1_3_3. There an assumed 
three times increase metalhcity of the main sequence progeni- 
tor of the white dwarf resulted in a ^^Ne mass fraction of 0.075. 
The 22Ne mass fraction was included in the post-processing by 
lowering the '^C fraction. 

With respect to the computational setup, model c3_3d dif- 
fers from the other three models. These were calculated on a 
static grid with uniform fine-resolved inner parts and exponen- 
tial grid spacing in the outer regions to capture at least parts of 
the expansion. Models b30, 2_2_2, and 1_3_3 were followed 
for 1.33 s, 1.5 s, and 2.0 s, respectively. The light curve deriva- 
tion assumes homologous expansion at the point where the hy- 
drodynamical explosion simulation ends. This is, however, not 
yet reached to a high accuracy at these times. Therefore, in 
model c3_3d a different approach was chosen (Ropke 2005). 
This simulation was carried out on a uniform computational 
grid co-expanding with the white dwarf. In this way it was pos- 
sible to follow the evolution time to 10 s after flame ignition. 

All simulations (except b30) were carried out on [256]'' 
grid cells to ensuring numerical convergence (Reinecke et al. 
2002a; Ropke 2005). To acconamodate for a reasonable num- 
ber of flame ignition kernels the resolution was increased to 
[756]3 ceUs for model b30. 

The global characteristics of the explosion models are sum- 
marized in Table 2. Detailed results from the explosion sim- 
ulations can be found in Travaglio et al. (2004), Ropke & 
HiUebrandt (2004, 2005b) and Ropke et al. (2005)^ AU mod- 
els produce rather weak explosions, lower than canonical 1 foe 
(= 10^' erg), and low amounts of ^^Ni. Although they he in 
the range of variability expected from observations, only model 
b30 comes close to the class of "Branch normals" (Branch et al. 
1993). 

^ Kozma et al. (2005) derived a nebular spectrum for model c3_3d. 
In their case this model was denoted c3_3d_256_10s. 
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The remapping of 3-D models to spherical 1-D geometries 
has been done using the same tracer particles that were used 
for the nucleosynthesis calculations. A grid in radial coordi- 
nates was constructed with uniform steps in radius, the outer 
radius being equal to the maximum radius of the 3-D model. 
The number of tracer particles (all of equal mass) then deter- 
mined the mass of each spherical shell constructed on the grid, 
and the chemical composition of the particles determines the 
composition of the shell. The square of the particle velocities 
were summed up to obtain the total kinetic energy of the shell. 
The motion of particles is not purely radial, however, the non- 
radial component is not high at the end of explosion simulation. 
The kinetic energy of transversal motion is of the order of 0.01 
of the radial motion (since model c3_3d was followed for a 
longer period the transversal motion is only ~ 3 x 10 '* of the 
radial motion). 

We assume that even if the non-radial component will be 
dissipated and goes to heat, the heat will not be radiated away 
because during the first several hours after the explosion all 
zones are optically thick. After some time this heat will be re- 
turned to the kinetic energy of the expansion. The kinetic en- 
ergy and the mass of each shell determines the velocity of the 
shell, which is ascribed to the outer radius of the mesh zone in 
the radiation hydro code. The radiation hydro simulation starts 
for all runs at time lO'* s after the ignition of the explosion. We 
assume that the total energy at the end of the flame simulation 
fitot (kinetic plus thermal plus gravitational) goes into at 
10^ s. We have renormalized the original Zikin from the tracer 
particles to Etot, and this value is given in the Table 2 as Ei_i„ 
initial. 

As an example, the chemical composition for the model 
2_2_2 is shown in Fig. 1. Here the mass fractions of the most 
abundant elements are presented as a function of the mass in- 
side each shell. This distribution of the composition is pre- 
served during the radiation hydro run for all elements (except 
for ^''Ni which decays successively to ^^Co and ^^Fe). A better 
presentation of the innermost and outermost layers is achieved 
when the composition is plotted as a function of the material 
velocity. This is shown for all models in Figs. 2-5. However, 
since not all of the initial models are in full homology, the ve- 
locity is not a good Lagrangian coordinate. One can notice in 
Fig. 2 that some outer shells have velocities lower than the ad- 
jacent inner shells. Moreover, the homology is never perfect 
during first weeks due to the heating from the ^^Ni — » ""^Co 
— » ^*Fe decay chain. Nevertheless, the plots in Figs. 2-5 do 
show the expected model distribution of chemical elements at 
the observable stages of supernova evolution with an accuracy 
of - 1%. 

Due to noise in the composition and density when averag- 
ing the 3-D model, the number of radial zones was restricted to 
50 (see, e.g. Fig. 1 and Fig. 6). No smoothing of the composi- 
tion or the density in the radial direction has been done during 
the 3-D to 1-D remapping. Note the value of any quantity in a 
spherical shell is just a mean value over An for the same shell 
in the 3-D model. 

Table 2 presents the parameters that are the most important 
for light curve modeling. These include: (1) the ^^Ni mass and 
(2) the initial and final kinetic energy. Note that the asymptotic 



Table 2. Mass of ^''Ni in the models and their energetics 



model 


M(^'^Ni)[Mo] 


iikin. foe, initial 


Ekin, foe, asymptotic 


1_3_3 


0.24 


0.357 


0.365 


2_2_2 


0.31 


0.441 


0.453 


c3_3d 


0.28 


0.431 


0.441 


b30 


0.42 


0.663 


0.679 




0.5 1 



Fig. 1. Mass fractions of the most abundant elements produced in 
model 2_2_2 as a function of Lagrangian mass. Note 'Fe' includes 
iron-peak elements together with ^^Ni. 

kinetic energy is somewhat higher than the initial kinetic en- 
ergy due to the addition of energy from the radioactive decay 
of ^^Ni. This small effect is discussed below (see Sect. 6.2). 

4. Predicted observables and their diversity 

Observations of SNe la indicate that a valid explosion 
model should release ~10^' erg of energy and synthesize 
~ 0.4 . . . 0.7 Mq of ^'^Ni. However, there exists a large di- 
versity in the observations ranging from the class of sub- 
luminous SNe la (SN 1991bg) to the bright events (SN 1991T). 
The deflagration models use in this study release 4.4 . . .6.3 x 
10^*^ erg of asymptotic kinetic energy into the ejecta and pro- 
duce 0.24 . . . 0.42 Mq of 5^Ni (see Table 2). Thus they fall into 
the lower range of observational expectations, but as of yet, do 
not account for the more luminous events. 

The synthetic light curves are sensitive to the energy re- 
lease, the ^''Ni production and the distribution of elements in 
the ejecta. Fig. 7 displays the computed B- and V-band light 
curves for all four of the explosion models (see Table 1). 
During maximum light these two passbands carry the main 
contribution to the luminosity. It is clear that in spite of modest 
values of explosion energy and ^^Ni mass, these deflagration 
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Fig. 2. Mass fractions of most abundant elements for the model 1_3_3 
as a function of the material velocity. 




Fig. 3. Mass fractions of the most abundant elements produced in 
model 2_2_2 as a function of material velocity. Note 'Fe' includes 
iron-peak elements together with ^'Ni. 



models quite nicely reproduce the observed absolute peak B- 
and V-band magnitudes and decline rates of observed weak and 
normal SN la. Below we compare the computed light curves 
with some well-observed nearby SNe la. 




Fig. 4. Mass fractions of most abundant elements for the model c3_3d 
as a function of the material velocity. Note 'Fe' includes iron-peak 
elements together with ""^Ni. 




Fig. 5. Mass fractions of most abundant elements for the model b30 as 
a function of the material velocity. 



4.1. Broad-band photometry 

Detailed comparisons with observations of SN 1992 A 
(Suntzeff, 1996) and SN 1994D (Richmond et al., 1995;Meikle 
et al., 1996; Smith et al., 2000) are presented in Figs. 8- 
11. These two SN la were selected because they are both 
well observed and the ^^Ni masses derived for each of them 
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Fig. 6. Initial (/ = 10" sec, black solid) and final (t = 90 days, dotted) 
density profiles as a function of material velocity scaled to the same 
maximum for c3_3d. 

(Stritzinger, 2005) are similar to the values of the models.^ The 
observational template light curves are obtained as described 
in Contardo et al. (2000) and Stritzinger & Leibundgut (2005). 
The filter transmission functions used to derive the light curves 
from the modeled spectrum are from Stritzinger et al. (2005). 
The zero time on the plots is for the time of B-band maximum 
(both for theory and observations). 

In Fig. 8 we compare synthetic light curves derived 
from model 1_3_3 (Ropke et al. 2005) with observations of 
SN 1994D. The distance to SN 1994D is still controversial. 
We have assumed the distance modulus 30.4 from (Drenkhahn 
& Richtler 1999) and the total reddening assumed is E(B- 
V)tot=0.022 

The model produced only 0.24 of ""^Ni, that is why 
the modeled UBV-bwd maxima are low compared with 
SN 1994D. The decline in B-band light curve for 15 days af- 
ter the maximum light is also slow for a low luminosity SN la. 
The near infrared maxima are of the correct order, although 
they lack a convincing secondary maxima. 

In Fig. 9 we compare computed light curves derived 
from model 2_2_2 (Ropke et al. 2005) with observations 
of SN 1992A. The distance modulus is taken to be 31.35 
and the observed light curves are de-reddened assuming E(B- 
V)tot=0.017. Richtler et al. (2000) have shown that a num- 
ber of different distance determinations of the Fornax cluster 
of galaxies (where SN 1992 A is located) agree well with a 
distance modulus of 31.35+0.04 mag (18.6+0.3Mpc). Madore 
et al. (1998) have published a Cepheid distance to the cluster 
giving a distance modulus of 31.35 + 0.07. 

^ The light curve templates of SN 1992A have less coverage in 
the R- and /-bands because there are fewer premaximum photomet- 
ric points compared to in the optical passbands. 



model 


tnvAB), days 


^max 


Am,5(fi) 


1_3_3 


17.58 


-18.00 


1.52 


2_2_2 


18.59 


-18.28 


1.40 


c3_3d 


16.96 


-18.26 


1.46 


b30 


18.18 


-18.73 


1.43 



The model produced 0.3 M© of Ni. This is not a sufficient 
amount to account for the peak B- and V-band magnitudes. 
However, the peak value in the f/-band is almost reached, al- 
though somewhat earlier than what is required by the observa- 
tions. The decline in the B- and V-bands and the shapes of I- 
and /?-band light curves are similar to those for model 1_3_3. 

Fig. 10 displays the theoretical light curves for model 
c3_3d as compared with observations of SN 1994D. The ''''Ni 
mass is rather low (0.28 Mq) and the general behavior of the 
fluxes is similar to the previous cases. 

The best agreement with observations in the B- and V- 
bands is achieved with model b30. Fig. 11 shows the syn- 
thetic light curves derived from model b30 with observations 
of SN 1992 A. Near maximum light there is a very good (al- 
most perfect) agreement in the B- and V-bands as well as in 
the decline rate 20 days past maximum. This is an important 
result for cosmological applications, since the decline rate in 
those passbands (which are the main contributors to luminos- 
ity at this epoch) is used for calibration of the absolute fluxes. 
The flux in the t/-band peaks early while at the same time over- 
shoots the observed peak luminosity, and then declines too fast. 
This behavior may be explained by the significant amount of 
mixing of ^^Ni in model b30 (see Fig. 5). With a large frac- 
tion of mixed ^^Ni, one would expect an excess of U -band flux 
prior to maximum followed by a rapid evolution. Less mixing 
of ^^Ni may help to put the model in a better agreement with 
observations. Note, however, the maximum U -band luminosity 
is on the correct order for SN 1994D. 

Our experiments with other models show that a more lay- 
ered structure of chemical composition may also help to pro- 
duce secondary maxima in the R- and /-bands, which are not 
reproduced by the strongly mixed model in Figs. 8-11. See 
Fig. 3. in Blinnikov & Sorokina (2004) for the abundance pro- 
file of the classical W7 model (Nomoto et al. 1984). 

Some computed light curves parameters are summarized in 
Table 3. If we compare our results for Amig(B) with the Phillips 
relation (Phillips et al. 1999), we find that model b30 fits quite 
well, while other lower ^''Ni mass models are somewhat under- 
luminous for their values of Ami5(B). In this comparison, one 
should take into account that we have used Ho - 72 km s"' 
Mpc while Phillips et al. (1999) used Ho = 65 km s ' 
Mpc"'. 

4.2. Bolometric light curves 

In Fig. 12 we compare the computed bolometric light curve for 
the model b30 with several of the Ultra Violet Optical InfRared 
(UVOIR) light curves presented in Stritzinger (2005). Note that 
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Fig. 7. B- and V-band hght curves for all 4 models. Left: zero time is 
the moment of explosion. Right: zero time is the maximum light in the 
given band. 



Fig. 9. UBVRI light curves for 2_2_2 compared with observed light 
curves from SN 1992A (data points). Zero time is S-band maximum. 




the theoretical bolometric light curve (dotted line) is computed 
in the wavelength range from extreme ultraviolet to far in- 
frared. The templates UVOIR light curves are computed in the 
range from f/-band (3100 A) to the near infrared (10,000 A). 
Therefore if we take a simple cut from the multi-group spec- 
trum and sum the flux we obtain a UVOIR light curve (solid 
line in Fig. 12) that can be directly compared to the template 
UVOIR light curves. Although the luminosity at maximum and 



the rise-time are in good agreement with SN 1992 A the decline 
rate is somewhat slower than observed. This is mostly due to 
the slow decline of the near infrared flux in the theoretical ra- 
diative model. 

In Fig. 13 we compare the bolometric light derived from 
model c3_3d to the template UVOIR light curves. Compared 
to the previous model, c3_3d has a lower explosion energy and 
produces a smaller amount of ^^Ni (see Table 2). As a results 
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Fig. 11. Synthetic light curves derived from model b30 (solid curves) 
(Travaglio et al. 2004) compared with observed light curves from 
SN 1992A (data points). Zero time is with respect to S-band maxi- 
mum. 



Fig. 12. Total bolometric luminosity (dotted) and UVOIR luminosity 
(solid) for the model b30. Observations, shown for three representative 
SNe la, should be compared only with the solid line. Zero time is the 
epoch of maximum luminosity for each light ciu've. 



we see that the total luminosity is less than what is determined 
from model b30. 

In the next section we compare the UVOIR luminosity 
for our four models to the luminosities derived from the tem- 
plate light curves by the method presented in Stiitzinger & 
Leibundgut (2005). 

4.3. Photospheric velocity 

Our radiation hydro models allow us to evaluate the velocity 
of material at the level of the photosphere. We define this level 
as a layer with an optical depth of 2/3 in the BV-bands. The 
results for the velocity at this level are shown in Fig. 15. We 
find that the photospheric velocity is near ~ 10^ km s ' and 
changes weakly in the models, while observations sometimes 
show larger velocities before maximum followed by a faster 
evolution. For example, several days before B-band maximum 
the Si II /16355 absorption feature in SN 1992 A indicated ve- 
locities on the order of 13 x 10^ km s"' (Kirshner et al. 1993). 
Moreover Patat et al. (1996) obtain for SN 1994D 14 x 10^ 
km s"' 10 days before B-band maximum. Similar values are 
found by Branch et al. (2005), see Fig. 15. In both SNe la 
the velocity fell quickly near maximum light and stayed -lO"* 
km s"' around day 20 after the maximum. Note, that in Fig. 15 
the time is given from the explosion epoch, and epochs of B- 
band maximum are presented in Table 3. 

It should be noted that Blondin et al. (2005) recently pre- 
sented a detailed analysis of line profiles for a large sample 
of SNe la. Their study indicates that there are many nuances 
at play that make the procedure of using line profiles as a mea- 
sure of photospheric velocities a difficult task. In addition some 
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Fig. 13. Bolometric luminosity (dotted) and UVOIR luminosity (solid) 
for the model c3_3d. 

SNe la exhibit a very slow evolution of velocity measured by 
Si II /16355 absorption near ~ 10"* km s ' , see Fig. 6 in (Blondin 
et al., 2005). 

This may be interpreted as the signature of lower bound of 
Si distribution (e.g.. Branch et al., 2005), not as the signature 
of the slow evolution of the photospheric velocity. To illustrate 
this we plot in Fig. 15 the behaviour of the velocity on the pho- 
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Fig. 14. Deposited gamma-ray (solid line) and bolometric luminosity 
(dashes) for the model b30. Zero time is the explosion epoch. 
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Fig. 15. Photospheric velocity for c3_3d (solid), 2_2_2 (dashed) and 
b30 (dotted) defined as /?t(sv)=2/3/?- Circles show the velocity at 

'Rosseland' photosphere for b30. Crosses are the data for SN 1994D 
from Branch et al. 2005. Zero time is the explosion epoch. 



ent level of excitation caused by the presence of Ni56 in outer 
layers. This may be enough to support optical depth in bluer 
wavelengths, but not sufficient for a 'true' photosphere. 

The observation of fast features in the spectra of some 
SNe la demonstrates the need for further development of the 
deflagration models. Our procedure of averaging of the original 
3-D models over solid angles somewhat suppresses the fastest 
motions. As a result the outermost mesh zone in model b30 
contains a mass of 8.5 x 10 '^M© with an asymptotic velocity 
12.98 X 10^ km s~\ while some of the tracer particles of the 
3-D model had 13.81 x 10^ km s'^ 



5. Deriving ^^Ni masses from SNe la ligiit curves 

Following "Arnett's rule" (Arnett 1979, 1982) one can derive 
the ^^Ni mass from the peak luminosity of a SN la. Arnett's 
rule simply states that at the epoch of maximum light the peak 
luminosity is equal to the rate of gamma-ray deposition inside 
the ejecta. The derivation of this statement is based on many 
simplifying assumptions, yet it is well satisfied in the current 
set of models (see Fig. 14). 

To derive a ^^Ni mass from the peak luminosity one also 
has to know the fraction of the total radioactive power that is 
actually deposited, and the fraction that escapes. An empiri- 
cal procedure for finding the ^^Ni mass has been developed by 
Stritzinger & Leibundgut (2005). Wifli a UVOIR light curve 
and the simple relation L,nax = 2 x 10"^^MNi/MO erg s"', (and 
10% correction taking into account the difference of UVOIR 
and bolometric luminosity) they were able to make estimates 
of the ^''Ni mass for a large number of SNe la. 

The accuracy of this procedure may be tested with our 
synthetic light curves. Fig. 16 contains the comparison of the 
UVOIR light curve derived from theoretical UBVRI fluxes ap- 
plying the procedure of Stritzinger & Leibundgut with the re- 
sults for model b30. 

Table 4 shows that the derived mass of ^''Ni is systemat- 
ically lower than the actual values in the models. Fig. 17 ex- 
plains the reason for this: the relation Lmax = 2 x lO'^^Mm/MO 
erg s"' is good for true (total) bolometric luminosity which is 
systematically higher than the UVOIR luminosity and the 10% 
correction is a bit small to fully account for that. 

The model b30 produced 0.42 Mq of ^^Ni and a compa- 
rable value of 0.40 M© of ^^Ni was found for SN 1992A by 
Stritzinger (2005). The agreement is better, than suggested by 
numbers in the Table 4, because the observed UVOIR luminos- 
ity at maximum fight is somewhat higher than in the model b30, 
since our models have IR fluxes too low at UVOIR maximum 
light (see Figs. Hand 12). 



tosphere determined by the Rosseland mean opacity for b30. 
While the the ejecta cool down after the maximum light the 
Rosseland opacity is determined by red wavelengths where the 
matter is more transparent than in blue. Here the photospheric 
velocity goes down monotonically which is usually observed 
in spectroscopic studies. The non-monotonic behaviour of the 
photospheric velocity in Fig. 15 may be explained by difl'er- 



6. Discussion of tlie resuits 

This study does not produce theoretical light curves that per- 
fectly match all the observations. However, the synthetic mod- 
eled UBV light curves do show good agreement with the ob- 
served ones for modestly luminous SNe la, although the evolu- 
tion of the luminosity is somewhat slower (particularly in the 
bolometric fight curve) than observed. In this section we dis- 
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Table 4. Values of ^''Ni mass derived from the theoretical light curves 
by the method from Stritzinger & Leibundgut (2005) . 



I, days 

Fig. 16. Computed UVOIR light curve for the model b30 (circles) 
compared with the UVOIR light curve derived using the procedure of 
Stritzinger & Leibundgut (2005) (sohd line) with the modeled UBVRI 
light curves. The full bolometric luminosity is given by the dotted line. 
Zero time is the explosion epoch. 
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Fig. 17. Peak values of UVOIR (stars) and bolometric (circles) lumi- 
nosity versus ^''Ni mass in logarithmic scale for our four models. Solid 
line is the relation from (Stritzinger & Leibundgut, 2005). Dashed ar- 
row points to the correct value of ''^Ni mass for h30. Thin solid line 
with arrow points to the lower ''^Ni mass found from UVOIR lumi- 
nosity and thick solid arrow is 10% correction used (according to 
Stritzinger & Leibundgut, 2005) to obtain the number given in the 
Table 4. 



model 


log L„,„, erg s ' 


M(5*Ni)[M0], 


M(5^Ni)[Mo], 






actual 


derived 


1_3_3 


42.579 


0.24 


0.21 


2_2_2 


42.681 


0.31 


0.26 


c3_3d 


42.636 


0.28 


0.23 


b30 


42.821 


0.42 


0.37 



cuss the approximations made in the current study and their 
probable influence on the results. 

6.1. Shortcomings of the current modeling 

Partially, the remaining differences between synthetic and ob- 
served light curves can be attributed to the the explosion mod- 
els. All centrally ignited models possess an artificial flame ig- 
nition configuration that is known to give rise to weak ex- 
plosions. In addition in the four explosion models used, ther- 
monuclear burning was not followed below densities of 10^ 
g cm"^. In a more realistic explosion model burning should 
continue to far lower densities. This has two effects that are 
likely to improve the agreement between the modeled and ob- 
served light curves. First, at lower densities more intermediate 
mass elements (IMEs) (e.g. Si, S, and Ca) will be produced 
which are clearly under-abundant in the cuiTent set of models. 
Second, burning to lower densities releases more nuclear en- 
ergy resulting in higher expansion velocities. Currently work is 
in progress to analyze the effects of extending the burning to 
lower fuel densities (see Ropke & Hillebrandt, 2005a). 

It is clear that the radiative transfer models need consider- 
able improvement at certain points. While up to 20 days after 
the maximum light we find good agreement of the best model 
with observations in the B- and V- bands, it is clear that at later 
epochs large deviations are present. During these epochs when 
the SN la ejecta enter the nebular phase, the LTE assumption is 
clearly not valid. In addition, the procedure of averaging a 3-D 
model for 1-D transport can introduce an error which can only 
be corrected in a full 3-D radiative transfer simulations. 

Another factor that has a dramatic influence on the modeled 
light curves is the treatment of the opacity. Direct numerical 
summation of individual lines, as done in atmospheric codes 
like PHOENIX (Baron et al. 1996b) may be dangerous in deep in- 
teriors of SN ejecta, since the difference of the source function 
and mean intensity is very small where the lines are strong. The 
subtraction of the nearly equal numbers always leads to loss of 
accuracy, and the result significantly depends on the order of 
summation of lines: one has to start from the weakest ones (cf. 
Sorokina & Blinnikov 2002). Opacity sampling (OS) or opac- 
ity distribution function (ODF) and ^-distribution procedures 
used in static atmospheres must be modified for applications in 
supernova envelopes (Wehrse et al. 2002). However, the recipes 
for the statistical treatment of lines were developed by Wehrse 
et al. (2002) only for the flux equation, and at present it is 
unclear how to generalize them for treatment with the energy 
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equation, e.g. Eq. (3). A detenninistic approach for an average 
expansion opacity to be used in the energy equation was pre- 
viously developed (Sorokina & Blinnikov 2002), but it is more 
computationally demanding compared to the simple formulae 
(Eastman & Pinto 1993) employed in this study. 

For some runs in the current set of models we have applied 
the new recipes by Sorokina & Blinnikov (2002) and found 
for them a small difference in comparison with the prescrip- 
tion by Eastman & Pinto (1993). This is probably due to strong 
mixing of the models. More layered models produce clear sec- 
ondary maxima in the R- and /- bands, and the use of 'new' 
opacities from Sorokina & Blinnikov (2002) makes them more 
pronounced. 

Recent papers by Branch et al. (2005); Stehle et al. (2005) 

used spectroscopic analysis in different approximations and 
found that the composition structure in SNe la is stratified but 
has some degree of mixing. The latter is not so strong as the 
mixing in the current models. This is another hint for need 
of producing more chemical stratification in the hydrodynamic 
models. However, one should first build detailed NLTE spec- 
troscopic models for the current models. The line spectra are 
sensitive not only to the distribution of chemical elements but 
also to the conditions of excitation which may change in radius 
in more pronounced way than the chemical composition. 

The absolute magnitudes for SNe la in the infrared J, H and 
A'-bands is very uniform even for pecuUar events (Krisciunas et 
al. 2001, 2004). This is a rather puzzling result given the variety 
of morphologies displayed in the R- and /- bands. The opacity 
physics in our radiation hydro model is presently not complete 
in the 7, H and /T-bands and it will be necessary to address this 
problem in future work. 

6.2. Hydrodynamic effects in SN la 

Progenitors of SNe la are believed to be compact degenerate 
stars. The explosion develops on a time-scale of a second, and 
homologous expansion should be a good approximation after 
a few seconds. This approximation is usually exploited in the 
radiative transfer codes that neglect hydrodynamics (Eastman 
& Pinto 1993; Lucy 2005). However, Pinto & Eastman (2000a) 
pointed out that the energy released in the ^^Ni decay can influ- 
ence the dynamics of the expansion. The ^^Ni decay energy is 
3 X lO'^erg g"' which is equivalent to a speed of 2.5 x 10^ km 
s if transformed to the kinetic energy of a gram of material. 
Pinto & Eastman (2000a) state "Since the observed expansion 
velocity of SNe la is in excess of 10"* km s"', we expect that 
this additional source of energy will have a modest, but perhaps 
not completely negligible, effect upon the velocity structure." 

In reality the heat released by the ^''Ni decay will not all 
go into the expansion of the SN la. If the majority of ^^Ni is 
located in the central regions of the ejecta then the main effect 
is an increase in the entropy and local pressure (both quanti- 
ties are dominated by photons in the ejecta for the first several 
week). The weak overpressure will lead to a small decrease 
in density at the location of the 'nickel bubble', as well as to 
some acceleration of matter outside the bubble. It is often said 
that the expansion of the ejecta is supersonic and pressure can- 



not change the velocity of the matter, but one should remember 
that in the vicinity of each material point we have a 'Hubble' 
flow so differential velocities are in fact subsonic in a finite vol- 
ume around each point. 

In the case of a more uniform distribution of ^*Ni, like in the 
models under consideration, we should not expect the forma- 
tion of a nickel bubble, but the general effect upon the velocity 
should be anticipated. 

In order to compute the hydrodynamics correctly to find 
small deviations from homologous expansion we should not 
assume that energy conservation can be simplified to the con- 
dition of thermal equilibrium, and that the gas terms dU + pdV 

nOO 

can be neglected. This simplification, when the term 4n L {rjy— 
ayJv)dv in the co-moving frame is equal to the rate of radioac- 
tive energy deposition s is always done in models that assume 
homologous expansion (Arnett 1980; Pinto & Eastman 2000a; 
Lucy 2005). 

We do not make simplifying approximations such as ne- 
glecting the gas terms dU + pdV, nor do we assume homolo- 
gous expansion. Eq. (3) automatically accounts for any devia- 
tions from thermal equilibrium. 

The effect is modest, but as expected it may result in a 
~ 1 0% difference in velocity. This effect is evident for the c3_3d 
model (see Fig. 6). The solid black line is the initial model 
scaled to our result at 90 days since explosion. We see that the 
density profile has changed due to the ^^Ni and ^^Co decays; 
it clearly deviates from homology. The growth of the velocity 
in outer layers is also visible. The effect upon the total kinetic 
energy is smaller (see Table 2). 

All of the models considered here have strongly mixed 
^^Ni. Other models, with modest mixing do not show appre- 
ciable increase of the velocity in outer layers, but they demon- 
strate that the nickel bubble, i.e. the depression of density in 
^*Ni-rich layers, continues to grow during the coasting stage. 
The change of the density is important for the deposition of 
gamma-ray energy, and this is reflected in the light curve. 

6.3. Secondary maxima in R, I. 

Most SNe la (except for subluminous SN 1991bg-like ones) 
display pronounced secondary maxima in their near-infrared 
and infrared light curves. Observations show a large variety of 
observed infrared maxima. According to Nobili et al. (2005), 
the light curve in /-band can peak before, as well as after B^ax 
(between -3 days and -i-3 days with respect to Bmax). On av- 
erage, the secondary maxima occurs 23.6 after the B^ax with 
a dispersion of cr = 4.4 days. However, there are strong ar- 
guments that this behavior is not random. Nobili et al. (2005) 
found a correlation between the time of the secondary peak and 
the B-band stretch parameter. A luminosity-width relation was 
also found between the peak /-band magnitude and the B- and 
/-band stretch parameter. Contardo et al. (2000) found some 
similar correlations between the strength of the secondary max- 
imum and times of derivatives of the light curve. 

This is a very interesting and long-standing problem for the 
theory of light curves. Although Eastman (1997) previously 
demonstrated that some models do produce a secondary in- 



S.I. Blinnikov et al.: Synthetic light curves of SNe la 



13 



frared maxima (in that case it was a sub-Chandrasekhar model), 
the theory is still far from providing robust predictions concern- 
ing the shape of /-band light curve and the correlations pointed 
out in Nobili et al. (2005). 

For example, Hoflich (1995) found a weak secondary max- 
ima for W7, and rather noisy /-band light curves for other mod- 
els. Pinto & Eastman (2000b) also found rather weak secondary 
maxima, depending strongly on the atomic line database, nev- 
ertheless Eastman (1997) and Pinto & Eastman (2000b) have 
presented a physical argument that attempts to explain the ob- 
served secondary maxima. According to them, aflcr ~ 20 days 
past maximum hght, when the monochromatic opacity in the 
near infrared becomes small due to recombination of higher 
ions, there exists a large amounts of singly ionized emitters 
(e.g. Ca II, Fe II, and Co II). Through fluorescence these lines 
cause an enhanced 'leakage' of previously trapped photons, 
thus giving rise to the secondary maxima. 

For this mechanism to produce a secondary maxima the 
SN la model must require a large mass fraction of Ca, Fe, and 
Co located in central layers that are emitting photons near 20 
days past maximum light. Too strong mixing of iron-peak ele- 
ments reduces their amount in the center. With the MPA models 
the total mass of Ca is relatively small. Model b30 has 0.0037 
M0 of Ca, while DD4 contains 0.0375 M© (Pinto & Eastman 
2000b) and W7 contains 0.025 M©; the latter two containing 
an order of magnitude more ! 

As previously mentioned, less pronounced mixing may al- 
low model b30 to produce R- and /-band Ught curves that agree 
with the observation more closely. With less mixing of Fe- 
group elements and IMEs like Ca, our models may produce R- 
and /- band secondary maxima. This may also lead to a delay 
in the f/-band rise-time. 

What should be changed in the deflagration models to de- 
crease the degree of the mixing, but at the same time to increase 
the amount of ^^Ni? We believe that the only feasible way is to 
lower the ignition density considerably. This will have two ef- 
fects. First, electron captures become unimportant, increasing 
the ^^Ni at the expense of the other NSE nuclei. Secondly, at 
low densities the '^C-^^C burning rate is low, too, and the heat 
diffusion time gets longer, thus reducing the laminar burning 
speed and giving the star more time to expand before rapid 
turbulent burning sets in. This should in fact also reduce the 
amount of mixing. These conjectures have to be tested in fu- 
ture detailed numerical modeUng. 

7. Conclusions 

How the diversity of local SNe la can aff'ect their use as distance 
indicators and the possibility of systematic trends (as a function 
of redshift) in their observed properties cannot be addressed 
without detailed radiative transfer modeling. 

The set of hydrodynamical models used in this studied syn- 
thesized modest amounts of ^^Ni ~0.3 - 0.4 Mq and displayed 
relatively low explosion energies (lower than the 'standard' 1 .5 
foe expected for a 'typical' SN la). However, our results show 
that Chandrasekhar mass models that burn by a pure deflagra- 
tion do produce both UBV light curves and photospheric ex- 
pansion velocities that match well with observed weak to nor- 



mal SNe la. As the majority of flux during maximum Ught is 
emitted within these passbands this is an encouraging result. 
Moreover, it is these passbands that are of great importance for 
cosmological appUcations of SNe la. It is clear that the mod- 
els require some improvement in explaining the shapes of the 
near infrared hght curves and in explaining fast spectral fea- 
tures that are observed in many normal events. The bolomet- 
ric light curves calculated from our deflagration models evolve 
slightly slower than what is indicated from observations. These 
discrepancies hint at the necessity of producing faster moving 
ejecta, and somewhat less mixed chemical composition. The 
latest deflagration (to be discussed in a subsequent publication) 
models are promising in this aspect. 

From the point of view of radiative transfer modeling sev- 
eral issues must be dealt with including: (1) determining the 
sensitivity of the results to the completeness of the atomic line 
data base, (2) the robustness of the various approximations em- 
ployed such as the expansion opacity prescriptions, and (3) de- 
termining the importance of (time-dependent!) NLTE effects. 
These issues are currently being addressed and small steps are 
being made to bring us closer to a full 3-D transport model for 
SNe la. 

Acknowledgements. This work was supported in part by the European 
Research Training Network "The Physics of Type la Supernova 
Explosions" under contract HPRN-CT-2002-00303. SB and ES are 
supported by MPA and UCSC guest programs and partly by the 
Russian Foundation for Basic Research (projects 05-02-17480, 04- 
02-16793). Special thanks to Nick SuntzefF for access to unpublished 
SN 1992A photometry. We thank the referee for his thorough report. 

References 

Arnett, W. D. 1979, ApJ, 230, L37 
Amett, WD. 1980, ApJ, 237, 541 
Arnett, W. D. 1982, ApJ, 253, 785 

Baron, E., Hauschildt, P. H., & Mezzacappa, A. 1996a, 

MNRAS, 278, 763 
Baron, E., Hauschildt, P H., Nugent, P., & Branch, D. 1996b, 

MNRAS, 283, 297 
Benetti, S., et al. 2004, MNRAS, 348, 261 
Blinnikov, S. I., Eastman, R., Bartunov, O. S., Popolitov, V. A., 

& Woosley, S. E. 1998, ApJ, 496, 454 
BUnnikov, S., & Sorokina, E. 2004, Ap&SS, 290, 13 
Bhnnikov, S. I., & Sorokina, E. I. 2000, A&A, 356, L30 
Blondin, S., Dessart, L., Leibundgut, B., et al. 2005, astro- 

ph/05 10089 

Branch, D., Fisher, A., & Nugent, P 1993, AJ, 106, 2383 
Branch, D., Baron, E., HaU, N., Melakayil, M., & Parrent, J. 

2005, PASP, 117, 545, 
Chan, K. W, & Lingenfelter, R. E. 1993, ApJ, 405, 614. 
Colgate, S. A., Petschek, A. G., & Kriese, J. T. 1980a, in 

AIP Conference Proceedings, No. 63: Supernova Spectra, 

ed. R. Meyerott & G. H. GiUespie (New York: American 

Institute of Physics), 7. 
Colgate, S. A., Petschek, A. G., & Kriese, J. T. 1980b, ApJ, 

237, L81. 

Contardo, G., Leibundgut, B., & Vacca, W. D. 2000, A&A, 
359, 876 



14 



S.I. BUnnikov et al.: Synthetic light curves of SNe la 



Drenkhahn, G. & Richtler, T. 1999, A&A, 349, 877 
Eastman, R. G. In: Thermonuclear Supernovae. ed. by P. Ruis- 

Lapuente et al. (Kluwer Academic Pub., Dordrecht, 1997) p. 

571 

Eastman, R. G. & Pinto, P A. 1993, ApJ, 412, 731 
Evans, K. F. 1998, J. Atmospheric Sciences, 55, 429 
Gamezo, V. N., Khokhlov, A. M., Oran, E. S., ChtcheUcanova, 

A. Y., & Rosenberg, R. O. 2003, Science, 299, 77 
HiUebrandt, W. & Niemeyer, J. C. 2000, ARA&A, 38, 191 
Hillebrandt, W., Reinecke, M., & Niemeyer, J. C. 2000, in 

Proceedings of the XXXVth rencontres de Moriond: Energy 
densities in the universe, ed. R. Anzari, Y. Giraud-Heraud, & 
J. Tran Thanh Van (Th§ Gioi PubUshers), 187-195 
Hoflich, P 1995, ApJ, 443, 89 

Hoflich, P. 2002, Workshop on Stellar Atmosphere Modeling, 
8-12 April 2002 Tuebingen. Eds: I. Hubeny, D. Mihalas, 
K. Werner, astro-ph/0207103. 
Hoflich, R, Khokhlov, A., & MuUer, E. 1993, A&A, 270, 223 
Hoflich, P., Wheeler, J. C, & Thielemann, F. K. 1998, ApJ, 
495,617 

Hoflich, R, Nomoto, K., Umeda, H., Wheeler, J. C. 2000, ApJ, 
528, 590 

Jefi^ery, D.L. 1999, astro-ph/9907015. 
Kirshner, R. R, et al. 1993, ApJ, 415, 589 
Kozma, C, Fransson, C, Hillebrandt, W, et al. 2005, A&A, 
437, 983 

Krisciunas, K., Philhps, M.M., Stubbs, C. et al., 2001, AJ, 122, 
1616 

Krisciunas, K., Philhps, M. M. & SuntzefF, N. B., 2004, ApJ, 
602L, 81 

Lucy, L. B. 2005, A&A, 429, 19 

Madore, B. E, et al. 1998, Nature, 395, 47 

Mannucci, R, DeUa Valle, M., & Panagia, N. 2005, MNRAS, 
astro-ph/0510315 

Meikle, W. P S., Cumnoing R. J., Geballe T. R., et al. 1996, 
MNRAS, 281,263 

Mihie, R A., The, L.-S., & Leising, M. D. 1997, in 
Proc. The Fourth Compton Symposium, ed. C. D. Dermer, 
M. S. Strickman, & J. D. Kurfess (New York: American 
Institute of Physics Press), 1022, astro-ph/9707111. 

Niemeyer, J. C. & Hillebrandt, W. 1995, ApJ, 452, 769 

Nobih, S., et al. 2005, A&A, 437, 789 

Nomoto, K., Thielemann, F.-K., & Yokoi, K. 1984, ApJ, 286, 

644 

Nomoto, K., Uenishi, T., Kabayashi, C., Umeda, H., Ohkubo, 
T, Hachisu, I. & Kato, M., 2003, in From Twilight to 
Highlight: The Physics of Supernovae, W. Hillebrandt & 

B. Leibundgut eds., ESO/Springer series "ESO Astrophysics 
Symposia" (Berlin), p.ll5 (astro-ph/0308138) 

Patat, R, Benetti, S., Cappellaro, E., Danziger, I. J., della Valle, 
M., Mazzah, R A., & Turatto, M. 1996, MNRAS, 278, 1 1 1 
Pignata, G., et al. 2004, MNRAS, 355, 178 
Pinto, P A. & Eastman, R. G. 2000a, ApJ, 530, 744 
Pinto, R A. & Eastman, R. G. 2000b, ApJ, 530, 757 
Phillips, M.M. 1993, ApJ, 413, L105 

Phiflips, M. M., Lira, P, Suntzeff, N. B., Schommer, R. A. 

Hamuy, M. & Maza, J., 1999, AJ, 118, 1766 
Pskovskii, Yu.P. 1977, Soviet Astron. 21, 675 



Reinecke, M., Hillebrandt, W, & Niemeyer, J. C. 1999a, 

Astron. Astrophys., 347, 739 
Reinecke, M., HiUebrandt, W., & Niemeyer, J. C. 2002a, A&A, 

386, 936 

Reinecke, M., HiUebrandt, W, & Niemeyer, J. C. 2002, A&A, 
391, 1167 

Richmond, M. W, Treffers, R. R., FiUppenko, A. V, et al. 1995, 

AJ, 109,2121 

Richtler, T, Drenkhahn, G., Gomez, M., & Seggewiss, W. 
2000, From Extrasolar Planets to Cosmology: The VLT 
Opening Symposium, Proceedings of the ESO Symposium 
held at Antofagasta, ChUe, 1-4 March 1999. Edited by 
JacqueUne Bergeron and Alvio Renzini. BerUn: Springer- 
Verlag, 259 

Ruiz-Lapuente, P. 1997, in Proc. NATO ASI on Thermonuclear 

Supernovae, ed. P. Ruiz-Lapuente, R. Canal, & J. Isem 

(Dordrecht: Kluwer), 68 1 
Ruiz-Lapuente, R, & Spruit, H. C. 1998, ApJ, 500, 360 
Ropke, F K. 2005, A&A, 432, 969 
Ropke, R K. & Hillebrandt, W. 2004, A&A, 420, LI 
Ropke, F. K., Gieseler, M., Reinecke, M. Travagho, C. & 

HiUebrandt, W. 2005, A&A, submitted 
Ropke, R K. & Hillebrandt, W. 2005a, A&A, 429, L29 
Ropke, R K. & HiUebrandt, W. 2005b, A&A, 431, 635 
Ropke, F. K., Hillebrandt, W, Niemeyer, J. C. & Woosley, S. E. 

2005, A&A, submitted 
Smith, C, et al. (private communication) 
Sorokina, E.I., BUnnikov, S.I. 2002. Nuclear Astrophysics, 

11th Workshop at Ringberg Castle, Tegernsee, Germany, 

February 11-16, 2002, ed. by E. MiiUer, W. ffillebrandt, 57 
Sorokina, E. I., BUnnikov, S. I., & Bartunov, O. S. 2000, 

Astronomy Letters, 26, 67 
Sorokina, E. 1., BUnnikov, S. 1., Kosenko, D. 1., & Lundqvist, 

P. 2004, Astronomy Letters, 30, 737 
Stehle, M., et al. 2005, MNRAS, 360, 1231 
Stritzinger, M. 2005, Technische Universitat Miinchen 

Dissertation 

Stritzinger, M., & Leibundgut, B. 2005, A&A, 431, 423 
Stritzinger, M., Suntzeff, N. B., Hamuy, M., Chalhs, P., 
Demarco, R., Germany, L., & Soderberg, A. M. 2005, RASP, 

117, 810 

Suntzeff, N. B. 1996, in: lAU CoUoquium 145, Supernovae 
and Supernova Renmants, ed. R. McCray, & Z. Wang, 
Cambridge University Press, Cambridge, 41 

Swartz, D. A., Sutherland, R G., & Harkness, R. R 1995, ApJ, 
446, 766 

Travaglio, C, Hillebrandt, W, Reinecke, M., & Thielemann, 

R-K. 2004, A&A, 425, 1029 
Wehrse, R., Baschek, B., & von Waldenfels, W. 2002, A&A, 

390, 1141 



